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Phase diagram of vortex states of high-Tc superconductors with sparse and weak columnar defects 
is obtained by large-scale Monte Carlo simulations of the three-dimensional anisotropic, frustrated 
XY model. The Bragg-Bose glass phase characterized by hexagonal Bragg spots and the diverging 
tilt modulus is observed numerically for the first time at low density of columnar defects. As the 
density of defects increases, the melting temperature increases owing to "selected pinning" of flux 
Unes. When the density of defects further increases, the transition to the Bose glass phase occurs. 
The interstitial liquid region is observed between these two glass phases and the vortex liquid phase. 

PACS numbers: 74.25. Qt, 74.62. Dh, 74.25. Dw 



In vortex states of high-Tc superconductors, columnar 
defects introduced by heavy-ion irradiation are strong 
pinning centers of flux lines, and dramatically enhance 
the critical current. [Q Nelson and Vinokur Q analyzed 
these defects and derived a phase transition between the 
Bose glass (BG) and vortex liquid (VL) phases. They 
also showed that this system can be mapped to two- 
dimensional interacting bosons in a random potential and 
that this phase transition corresponds to a localization 
transition of interacting bosons. In the case of dense 
and strong columnar defects, every flux line is trapped 
by randomly-distributed columnar defects at low tem- 
peratures, and the BG phase naturally appears. Many 
numerical studies ^ ^, p[ have been made in this case. 

In the case of sparse and strong columnar defects, part 
of flux lines are not trapped by defects, and physical 
properties are more complicated. At low temperatures 
untrapped flux lines are also frozen owing to interactions 
with trapped flux lines, and the BG phase is stable. At 
intermediate temperatures the defects still trap some flux 
lines but other interstitial ones are freely moving around, 
and the resistivity becomes finite. This region is called 
as the interstitial liquid (IL). At high temperatures all 
flux lines are depinned and the VL phase appears. It is 
not clarified yet whether there exists a phase transition 
or merely a crossover between the IL and VL. A numer- 
ical simulation was attempted in this case, and the 
existence of the BG phase was confirmed. Two kinds of 
BG phases for dense or sparse columnar defects are di- 
vided by the Mott insulator phase at the matching field. 
Recently the transition between these two (strong and 
weak) BG phases was observed numerically. H] 

In pure systems, the triangular fiux-line lattice (FLL) 
is formed below the melting temperature Tm- Is such 
order completely destroyed by introducing infinitesimal 
columnar defects? In vortex states with point defects, 
the long-range order of vortex positions is destroyed by 



infinitesimal defects, but the Bragg glass phase [l^ 
characterized by the quasi long-range order of vortex po- 
sitions and Bragg spots of the structure factor is sta- 
ble up to a considerable density of point defects. In the 
present Letter, we numerically show that a similar phase 
with Bragg spots associated with hexagonal symmetry is 
also observed in the case of sparse and weak columnar 
defects. Since the tilt modulus shows a diverging behav- 
ior in this phase, it corresponds to the Bragg-Bose glass 
(BBG) phase. Jill When the density of columnar defects 
increases, a phase transition to the BG phase takes place. 
We also find that the IL region exists above the BBG or 
BG phases. These results are summarized in Fig. |l|. 

In this study we perform large-scale Monte Carlo simu- 
lations based on the three-dimensional anisotropic, frus- 
trated XY model on a simple cubic lattice na] : 
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where the magnetic field is applied along the c axis. The 
periodic boundary condition is applied in all directions. 
A columnar defect is defined as a set of four weak inter- 
actions Jij = (1 — e) J surrounding a plaquette in the ab 
plane. Such defects distribute randomly in the ab plane 
with probability p, and take the same position along the 
c axis. Other interactions in the ab plane are given by 
Jij = J. This modeling of columnar defects is a vari- 
ant of our previous study on point defects. JIJ] Here we 
concentrate on the model with the anisotropy F = 5, the 
density of fiux lines / = 1/25, the strength of defects 
e = 0.1, and the system size T^ = T^ = 50 and Lc = 80. 
CPU time is 1 - 6 x 10^ Monte Carlo steps (MCS) for 
equilibration, and 1 ~ 3 x 10^ MCS for measurements. 
Since such long simulation time is required for taking 
equilibrium data, the present results are based on one 
typical configuration of columnar defects for each p. 



2 



M-l 



-I — I — I — I — I — I — I — I — ^ — I — I — I — I — ^ 
/ 



oel Bose glass /interstitial/ 
liquid ^ 

/ vortex- 
/ liquid 



iBragg-Bose 
glass 




T [J/k] 

B 



FIG. 1: Phase diagram of vortex states of high-Tc supercon- 
ductors with sparse and weak columnar defects in the plane of 
number of columnar defects per flux line versus temperature, 
p// = 1 corresponds to the matching field. 



First, we fix the density of columnar defects for p = 
0.01 {p/f = 0.25) and vary the temperature. From now 
on, this density is expressed in terms of the number of 
defects per flux line p/f with / = 1/25 in order to clarify 
the physical situation. The specific heat C shows a sharp 
peak around T ^ 0.37J//cb (marked in Fig. ||(a)), and the 
maximum value of the structure factor drops at the same 
temperature (Fig. §(b)). These facts indicate that the 
melting transition occurs at Tm ~ 0.37J//cb. The struc- 
ture factor shows sharp hexagonal Bragg spots below 
as displayed in Fig. |2|(c). Coexistence of two phases is ob- 
served between T = 0.368 J//cb and 0.374 J//cb owing to 
finite-size effects. These findings suggest that the transi- 
tion at Tm is of first order. Similar results are obtained 
for p/f = 0.075, 0.125 and 0.30. 

Since the helicity modulus along the c axis Tc takes a 
finite value when flux lines are frozen along the c axis, the 
vanishing point of Tc corresponds to the boundary be- 
tween the IL and VL. As shown in Fig. |^(a), temperature 
dependence of this quantity changes around Tm, and it 
vanishes at Til ^ 0.392 J//cb, which is close to the broad 
peak of C. Simulations for p/ f = 0.075, 0.125, 0.50 and 
1.00 show similar temperature dependence of Tc except 
for monotonic increase of Til as p increases. Size depen- 
dence of Tm and Til is checked for p/f = 0.075 and 0.25 
using data with Lc = 54, 80 and 160, and found to be 
smaller than that of T^ of the pure system. 

As displayed in Fig. §(b), the ratio of entangled flux 
lines to total flux lines A^ent /Afflux rapidly increases from 
0.3 to 0.5 around T^ as in the pure system. [|l^ Grad- 
ual increase of A^ent /Afflux in the IL region indicates the 
process that untrapped flux lines are freed from trapped 
ones. A'ent /Afflux — 0.85 at Til is comparable to the value 
at the onset of the VL phase in the pure system. [|l5| 

Second, we vary the density of columnar defects for 
T = 0.300J//cb. As shown in Fig. ||, the maximum 
value of the structure factor drops at the BBG-BG phase 
boundary. Although anomalies of thermodynamic quan- 
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FIG. 2: Numerical results for p/ f — 0.25. (a) Temperature 
dependence of the specific heat (triangles) and the helicity 
modulus along the c axis (squares) . Marked is the sharp peak 
of the specific heat at Tm- (b) Temperature dependence of 
the maximum value of the structure factor (circles) and the 
ratio of entangled flux lines to total flux lines (diamonds). 
Coexistence of two phases is observed in the hatched region, 
(c) Structure factor at T = 0.36J/A;b. 



titles such as the internal energy are too weak to observe 
numerically in the present study, the change of the sym- 
metry of Bragg spots is quite apparent as depicted explic- 
itly in Fig. ||. The hexagonal symmetry of Bragg spots 
at p/f = 0.35 (BBG phase) is lost at p/f = 0.375 (BG 
phase) and a glassy ring-like pattern appears. Similar 
results are obtained for T = 0.330 J//cb and 0.360J//cb. 

Third, in order to evaluate the BG transition temper- 
ature Tg, the tilt modulus C44 is measured. The in- 
verse of this quantity is proportional to the zero-field 
susceptibility of the transverse magnetization for the 
transverse field H'^ , namely c^^ ~ 
The Hamiltonian under the transverse field is given by 
H = Ho — H^m^, where Ho stands for the Hamilto- 
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FIG. 3: Defect density dependence of the maximum value of 
the structure factor for T = 0.300J//cb. Structure factors at 
p/f = 0.35 and 0.375 are also displayed. 
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FIG. 5: Top view of averaged positions of flux lines with 
respect to 80 planes and 2000 independent snapshots at (a) 
p/f = 0.25 and T = 0.360 J/Zcb (BBG phase), (b) p/f = 0.25 
and T = 0.380 J/Zcb (IL region), and (c) p/f = 0.40 and 
T = 0.360 J//cb (BG phase). Columnar defects are painted 
in yellow (no flux lines) to red (fllled with flux lines), and 
other regions are painted in violet (no flux lines) to light blue 
(largest probability of flux lines). 
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FIG. 4: Temperature dependence of the tilt modulus deflned 
in eq. (||) for (a) p/f = 0.125 to 0.30 and for (b) p/f = 0.35 
to 1.00. 



nian without the transverse field given in eq. (|l|) , and the 
transverse magnetization is defined from the sum of 
the difference of positions of flux lines along the direction 
of at the top and bottom ab planes. Explicitly, C44 
is given, aside from the coefficient, by 0] 



C44 



(3) 



Temperature dependence of C44 for p/ f = 0.125 to 0.30 
(BBG-IL boundary) and for p/f = 0.35 to 1.00 (BG-IL 
boundary) is displayed in Figs, ^(a) and ^(b), respec- 



tively. In both cases, C44 shows a large jump (note that 
these figures are drawn in semi-log scale) at Tg. This 
jump is expected to be infinite in the thermodynamic 
limit. Since Tg coincides with for p/f< 0.30 (see 
Figs. and ^(a)), the diverging behavior of C44 is ob- 
served in the whole region of the BBG phase. 

Next, real- space images of flux lines are shown in Fig. 
^. Positions of flux lines are averaged in all ab planes and 
during 10^ MCS (measured for every 500 MCS). Details 
of these rainbow-color figures are explained in the cap- 
tion. At p/f = 0.25 and T = 0.360J/A:b (Fig. |(a)), all 
flux lines are localized and seem to form a triangular lat- 
tice, which characterizes the BBG phase. At p/ f = 0.25 
and T = 0.380J/A:b (Fig. ^h)), the overall order does 
not exist any more, though columnar defects still trap 
flux lines. These properties are consistent with those of 
the IL region. At p/f = 0.40 and T = 0.360J/A:b (Fig. 
^(c)), although all flux lines are localized, the FLL struc- 
ture seen in Fig. ||(a) is totally destroyed. Since most 
columnar defects are occupied by flux lines unless other 
occupied ones are nearby, many dislocations inevitably 
appear. These properties signal the BG phase. 

Finally, interesting points of the phase diagram given 
in Fig. |l| are discussed. The most important feature is the 
existence of the BBG phase, which has not been observed 
experimentally or numerically until present. Giamarchi 
and Le Doussal [|l6| argued that nature of the glass or- 
der of the three-dimensional vortex states with columnar 
defects is equivalent to that of the two-dimensional vor- 
tex states with point defects, and pointed out |ll| the 
possibility of the BBG phase. Stability of the quasi long- 
range order against topological dislocations was discussed 
by D. S. Fisher, and he showed that this order is de- 
stroyed in d = 2 but stable in d = 3. Our finding is 
consistent with these analytic studies. 

The melting temperature of the BBG phase in- 
creases as the density of columnar defects p increases 
up to p/f = 0.30. This behavior may look strange at 
a glance because the FLL order is enhanced by intro- 
ducing randomness. The key to understand this exotic 
behavior is the real-space image of fiux lines in the BBG 



4 



phase given in Fig. ||(a). There exists some unoccupied 
columnar defects (yellow dots) even though no occupied 
defects are nearby. This fact suggests that the elastic re- 
pulsion between flux lines is dominant in the BBG phase 
and makes the columnar defects far from stable positions 
of FLL remain unoccupied. Owing to this "selected pin- 
ning" mechanism, increase of p does not seriously affect 
energy cost for FLL distortion, but bring energy gain 
from pinning of flux lines by the columnar defects close 
to stable positions of FLL. In this case, increases as p 
increases. As temperature increases, the BBG melts into 
the IL and the number of occupied flux lines increases as 
shown in Fig. |(b). This behavior can also be explained 
naturally by the above mechanism. 

When p further increases, energy gain from pinning 
of flux lines by many columnar defects exceeds that from 
FLL formation, and the transition to the BG phase takes 
place. The position of the BBG-BG phase boundary thus 
depends on the strength of columnar defects e. When e 
increases, this phase boundary shifts toward smaller p. 
Actually, simulations for e = 0.5 p/ f = 0.25 results in 
the BG phase at low temperatures, where the BBG phase 
exists for e = 0.1 (Fig. ||). There might be a critical value 
ec above which no BBG phase exists for any p. 

The boundary between the BBG or BG and IL phases 
(on which C44 jumps) is not monotonic in the vicinity of 
the BBG-BG phase boundary. This result shows a sharp 
contrast to the cage-model analysis, [|l^ where the "melt- 
ing temperature" monotonically increases as the strength 
of columnar defects increases. In the cage model, the 
system is reduced to a single flux line in a mean-field 
potential, and the BBG and BG phases cannot be dis- 
tinguished any more. We consider that the monotonic 
melting line obtained from the cage model corresponds 
to the BG-IL phase boundary above p/ f = 0.40. When 
this smooth curve is extrapolated to p = 0, it seems to 
meet the melting temperature of the pure system, T^p^^®. 
The excess superconducting region below p/ f = 0.40 is 
due to a nontrivial interplay between columnar defects 
and FLL structure, which is the origin of the BBG phase. 



In the melting transition of FLL in pure systems, sud- 
den entanglement of flux lines and increase of dislocations 
in the ah plane occur simultaneously, and dominance of 
these two contributions has been discussed. Random 
columnar defects along the c axis suppress entanglement 
of flux lines while enhance dislocations in the ah plane. 
The relation > T^p^^® might suggest that proliferation 
of entanglement of flux lines rather than that of dislo- 
cations in the ah plane is the main driving force of the 
melting transition. 

Three-dimensional vortex systems with columnar de- 
fects can be mapped to the two-dimensional interacting 
bosons in a random potential in the ground state. ^ 
Then, the present result indicates that the quasi long- 
range order of interacting bosons in two dimensions is 
stable in the presence of sparse and weak randomness. 

In conclusion, the phase diagram of vortex states of 
high-Tc superconductors with sparse and weak columnar 
defects is obtained by large-scale Monte Carlo simula- 
tions as shown in Fig. The Bragg-Bose glass phase 
characterized by Bragg spots with hexagonal symmetry 
and the diverging tilt modulus is observed for the first 
time. The melting temperature to the interstitial liquid 
region increases as the density of columnar defects in- 
creases owing to a "selected pinning" mechanism. When 
the density of defects further increases, a transition to 
the Bose glass phase takes place. The interstitial liquid 
region lies between the Bragg-Bose glass or Bose glass 
phases and the vortex liquid phase. 
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